Kondo Cloud Mediated Long Range Entanglement After Local Quench in a Spin Chain 
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We show that, in the gapless Kondo Regime, a single local quench at one end of a Kondo spin chain induces a 
fast and long lived oscillatory dynamics. This quickly establishes a high quality entanglement between the spins 
at the opposite ends of the chain. This entanglement is mediated by the Kondo Cloud, attains a constant high 
value independent of the length for large chains, and shows thermal robustness. In contrast, when the Kondo 
cloud is absent, e.g. in the gapped dimer regime, only finite size end to end effects can create some entanglement 
on a much longer time-scale for rather short chains. By decoupling one end of the chain during the dynamics one 
can distinguish between this end-end effect which vanishes, and the global Kondo cloud mediated entanglement, 
which persists. This quench approach paves the way to detect the elusive Kondo cloud through the entanglement 
between two individual spins. Our results show that non-perturbative cooperative phenomena from condensed 
matter may be exploited for quantum information. 
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I. INTRODUCTION 



Entanglement in many -body spin systems is a topic of high 
interesti— where "quantum correlations" or long range en- 
tanglement between individual spins is notoriously difficult 
to achieve. Here we propose to generate such entanglement 
through a non-perturbative dynamical mechanism which re- 
quires only a minimal action on a spin chain, namely a sud- 
den quench of a single bond. We consider spin chains with a 
magnetic impurity 8 located at one end; under certain condi- 
tions, the ground state of the chain is such that the impurity is 
maximally entangled to the block of spins forming the Kondo 
Cloud (KC)-2&. We will show that, during the dynamics in- 
duced by the quench, it is the Kondo cloud which mediates 
entanglement between the two ending spins of the chain. The 
generated entanglement is measurable through two spin corre- 
lations and enables a direct detection of the presence of a KC 
within the system before quenching. This is important since, 
despite several manifestations, the direct observation of a KC 
is still an open challenge^. Besides being a quantum infor- 
mational tool for probing a condensed matter paradigm, the 
long range entanglement created by local quenching is sub- 
stantial and potentially useful for linking separated quantum 
register*^. 

Though entanglement in condensed matter systems is typi- 
cally very short ranged 2 , there have been a few other propos- 
als for long distance entanglement - however they come at a 
high price. For example, there are proposals exploiting weak 
couplings of two distant spins to a spin chain 12 ' 13 , which have 
limited thermal stability or very long time-scale of entangle- 
ment generation. Otherwise, a dynamics has to be induced 
by large-scale changes to the Hamiltonian of a systemic, as 
opposed to the minimal change required in our proposal. Ad- 
ditionally, the exploitation of a non-perturbative cooperative 
feature of a condensed matter system, such as a Kondo cloud, 
has not been addressed so far. 

Quantum spin chains allow to investigate non-perturbative 
phenomena arising from the presence of impurities 1 -. For in- 
stance, the spin chain Kondo mode l 6 ' 7 ] 16 arises when a mag- 
netic impurity^ is coupled to the end of an antiferromagnetic 



Ji — J2 spin- 1/2 chain, where Ji(Ja) is the (next) nearest 
neighbor coupling. The model supports a crossover from a 
gapless Kondo Regime (KR) for J 2 < Jf = 0.2412 J x to a 
gapped Dimerized Regime (DR) for J2 > Jf. It is called a 
Kondo model since the impurity spin forms, in the KR, an ef- 
fective singlet with a the spins inside the KC. Recently, the en- 
tanglement of the ground state of this system has been charac- 
terized by varied means 6 ^ 7 .. The Kondo spin chain is described 
by the Hamiltonian 
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where cr, = (o; 
i, Ji > (J a 



, a\ , of ) is the vector of Pauli operators at site 
> 0) is the anti ferromagnetic (next) nearest 
neighbor coupling constants (we always put J\ = 1), J' < 1 
is the impurity bond strength and N is the number of spins, 
which we assume to be even through this paper. Hi is the ini- 
tial Hamiltonian of our dynamical scheme. Tuning J2 enables 
one to shift the spin chain between the Kondo and DRs. 



II. LONG-RANGE DISTANCE-INDEPENDENT 
ENTANGLEMENT IN THE KONDO REGIME 

We consider the finite Kondo chain ([T]i in its ground state 
\GSi). We, then, pertinently quench the coupling at the oppo- 
site end of the impurity allowing for the dynamics to develop 
entanglement. We show that, in the KR, the entanglement be- 
tween the two ending spins oscillates between high peaks with 
a periodicity determined by J', while the dynamics is very fast 
(thereby decoherence hardly gets time to act) and is robust 
against thermal fluctuations. In the dimer phase, the dynam- 
ics is much slower, qualitatively different and, in finite chains, 
it generates some entanglement due to unavoidable end to end 
effects, which are drastically tamed if one "cuts off the im- 
purity from the chain during the dynamics. In the KR, cutting 
off the impurity has minimal effect on the final entanglement 
between the ending spins since, here, the process is mediated 
by the cloud. 
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Entanglement is expected to be very different&i2 as J2 
crosses J| since, for J 2 < Jf , there is a characteristic length 
- the so called Kondo screening length £ - determined only by 
J' through £ oc e"'^ where a is a constant. In the KR, for 
any given J' and N, £ determines the size of a domain whose 
spins are maximally entangled with the impurity spin sitting 
at the origin^. 

Initially, the system is assumed to be in the ground state 
\GSi) of Hi. A minimal quench modifies only the couplings 
of the Ath spin by the amount J' (same as J' in Eq.(Q]i) so 
that Hi is changed to 

Hp = J\J\0\.02 + Jl<J\-Oz + J\<jn-\-<tn + J2&n~2-&n) 
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Since \GSi) is not an eigenstate of Hp it will evolve as 
\tP(t)) = e~ lHFt \GSi). An entanglement E(N,t,J') be- 
tween the ending spins emerges as a result of the above evo- 
lution. To compute E(N,t, J'), we first obtain the reduced 
density matrix piN(t) = tr 1 j v |^(t))(V'(t)| of spins 1 and 
A by tracing out the remaining spins from the state \ijj(t)). 
Then, we evaluate E(N, t, J') in terms of a measure of entan- 
glement valid for arbitrary mixed states of two qubits called 
the concurrence-^. Entanglement takes its maximum E m at 
an optimal time t opt and an optimal coupling J' opt such that 
E m = E(N,t op t, J' pt )- As we shall see J' opt is not a per- 
turbation of J 1 and J 2- If, as expected from scaling in the 
KR2ii£, the dependence on N and t can be accounted for by 
a redefinition of J' (equivalently £), then t opt and J' t cannot 



be independent quantities 
deed that, in the KR 



Our numerical analysis shows in- 



t opt oc A and that J' opt yields £ = N-2; 



since £ ex e Q/v J one gets t opt oc A oc 

For our choice of J', J\ and J 2 the spin-chain dynamics is 
not analytically solvable and one has to resort to numerical 
simulations. Recent methods of many-body simulations al- 
low handling exponentially big Hilbert spaces with pertinent 
truncations. Here, for N > 20, we use the time-step targeting 
method, based on the DMRG algorithm introduced ini^. For 
A < 20, instead, we resort to exact diagonalization. 

We focus only on the first period of the entanglement evo- 
lution in both regimes, since both decoherence and numerical 
errors make it unwise to wait for longer times. Fig. [Tfa) shows 
that fast long-lived (non-decaying) periodic oscillations with 
a period of 2t opt characterize the time evolution in the KR and 
that the maximal entanglement is achieved when the impurity 
coupling J' equals the value J' opt associated to a KC of size 
£ = A — 2 (the KC generated by the impurity sitting on the 
left side touches the other side of the chain); in the DR the 
dynamics appears more dispersive and not oscillatory for any 
J'. In Fig. [TJb) we plot the maximum of entanglement, E m , 
induced by bond quenching as a function of the length A: 
though the entanglement decreases as A increases, its value, 
in the KR, stays rather high and becomes almost distance in- 
dependent for very long chains; furthermore, as A increases, 
the entanglement generated in the KR is significantly bigger 
than the one in the DR. Despite its lower value, achieving en- 
tanglement in the DR costs more time, as shown in Fig. [TJ C )- 




FIG. 1: (Color online) Comparing the Kondo (J2 = 0) and dimer 
(J2 = 0.42) regimes, a) Entanglement vs. time for N = 30 with 
J' = 0.19 for the KR (solid line) and J' = 0.44 for the DR (dashed 
line), b) E m vs. the length N. c) t op t vs. N. d) J' opt vs. length N. 



It is also clear from Fig. [TJc) that t opt increases by A lin- 
early. Finally, in the KR, J' opt slowly decreases as A increases 
while it stays essentially constant in the DR (Fig. |TJd)). This 

is commensurate with the expectation that t opt oc e a ^ J ° pt in 
the KR, while, in the DR, t op t and J' opt are two independent 

quantities. The plot in Fig. |2K shows the exponential depen- 

z 1/2 

dence of t opt on J opt 1 realized, in the KR, for long enough 
chains. 

How the dynamics creates - even for very long chains of 
size N- high entanglement oscillations between the ending 
spins of the chain in the KR? To understand this, one should 
recall that, in the KR, the impurity spin forms an effective sin- 
glet with all the spins inside the cloud- and that, only in this 
regime, one can always choose J' so that £ may be made com- 
parable with N; at variance, in the DR, the impurity in \GSi) 
picks out - no matter what the values of J' and N are- only 
an individual spin in the chain to form a singlet ( a valence 
bond) while the remaining spins form singlets (local dimers) 
with their nearest neighbors^. Thus, only in the KR, one may 
use the remarkable resource to select - for any N- an initial 
state \GSi) which is free from local excitations: in fact, when 
£ = A — 2 (J' = J ' opt ) there is only a single entity, namely 
the Ath spin, interacting with the impurity-cloud composite. 
As this situation can always be engineered by choosing, for 
any N, J' = J' opt , this explains the distance-independent en- 
tanglement in Fig. QJ. The proposed scenario provides an in- 
tuitive grasp on why, only in the KR, an optimal entanglement 
between the ending spins may emerge from quench dynamics 
as the result of the interplay between very few states. At vari- 
ance, in the DR, the energy released by quenching is dispersed 
over the variety of different quantum modes associated to the 
local dimers. 

To provide a more quantitative analysis, one may expand 
|?/>(i)) in terms of eigenvectors of Hp. By exact diagonal- 
ization (up to N = 20), one finds that, in the KR, only two 
eigenstates of Hp (the ground state \E\) and one excited state 
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I-E2}) predominantly contribute to the dynamics: 

\Eh) = a 1 |^-) 1Ar |0-) 6 + /3 1 (|OO) 1A r|0 oo ) b 

+ \ll) 1N \(l> n ) b -\ij+)^+) b ) 

\E 2 ) = a 2 \^) lN \<j)-) b - p 2 {\00)iN\4> m )b 

+ lll)^^ 11 )^-!^)!^)^). (3) 

In Eq. ([3]l the first and the last spin are projected onto the 
singlet (|V' _ )) and the triplets (|00), 1 11) and |?/> + )) while the 
states of all spins in the body of the chain have been specified 
by the index b. After a time t the state evolves as \ip(t)) = 
(£i|GS7)|£i) + e- lAEt (E 2 \GS I )\E 2 ) + where AE is 
the energy separation between the two levels. One defines t = 
t op t as the time for which the contribution of |V'~)iA r l ( / ) ~)b 1S 
most enhanced in \ip(t)} due to a constructive interference. 
The condition for the onset of constructive interference is 

|( J B 1 |G5 / )ft|«|(^ 2 |G5 J )/3 2 |, (4) 

so that terms other than \il)~)iN\4>~)b m IV-'W) do not con- 
tribute at t = t opt . When J' — > 1 (very small cloud) then 
\{E!\GSi)\ w 1 while \{E 2 \GSi)\ as the ground state is 
hardly changed on quench: thus, constructive interference be- 
tween \Ei) and \E 2 ) it is impossible. The condition of Eq. © 
cannot be satisfied also when J' — > 0, since one has now that 
/3j ~ and f3 2 w 1 (the end spins form a singlet and triplet 
with each other in \Ei) and |J5a) respectively 12 ). Thus, only 
for intermediate J' entanglement may peak. When J' > J opt 
(£ < A — 2)- and particularly for £ < AT/2- the entanglement 
between the ending spins is frustrated by the existence of lo- 
cal excitations whose number increases as the size of the cloud 
gets smaller. In addition, when J' < J opt (£ > N — 2), the 
KC overtakes the chain and the Nth spin is already included 
in the cloud and its tendency is to screen the original impurity 
as in \4>~)iN\4>~)b rather than to pair with it to form a spin 
one (as in the last three terms of \Ei)). This makes /3i quite 
small, and it becomes smaller as the cloud overtakes the chain 
and again the condition of Eq. cannot be fulfilled. Con- 
sequently, the optimal situation is realized when J' = J' t 
(£ = N — 2 ), i.e just before the cloud overtakes the chain. 
Thus, only in the KR, one can convert-for any A- the useless 
entanglement between the impurity spin and the KC into a us- 
able entanglement between the ending spins of the chain. The 
emerging long distance entanglement analyzed in this paper 
is, indeed, a genuine footprint of the presence of the KC in 
\GSi) (Fig. Ha)). 

In Fig. [2J5 we plot- for both regimes- the entanglement 
reached after waiting for a time interval of the order of 1/ J' opt . 
One notices that, for J 2 > Jf , the entanglement peak de- 
creases sensibly and goes to zero rather soon. The plot of the 
maximal entanglement vs. J' is given in Fig. |2j;: here one 
sees that, in the KR, the entanglement rises from zero already 
at very small values of J'. This is expected since, in the KR, to 
a small J' is associated a large cloud containing the impurity 
sitting on the left side. 

The essential role of the KC in the entanglement generation 
is further probed if one let evolve the ground state \GSi) with 
a doubly quenched Hamiltonian obtained from (f2]i by isolat- 
ing (i.e., putting J' = 0) the left hand side impurity while 




j' 



FIG. 2: (Color online) a) t opt vs. l/y/J' opt in the KR. b) E m at 
t oc 1/ Jopt vs. J2 for chains of different lengths, c) Entanglement 
attained at t opt vs. J' for Kondo and DRs for N = 20. d) E m vs. 
temperature after bond quenching (blue line) and induced by (static) 
weak coupling with the rest of the chain (red line) for a chain of 
N = 10. 



keeping fixed to J' opt the bond connected to the last spin (see 
Fig. [3})). This forbids the dynamical build up of that "por- 
tion" of the total entanglement which is only due to end to 
end effects. In Fig. [4] we have plotted the E m vs. A after 
double quenching in both regimes. Fig. [4] shows that entan- 
glement in the dimer phase collapses already when A > 12 
while it stays unexpectedly high - and almost independent on 
A - in the KR; the existing entanglement of the KC with the 
impurity- is dynamically swapped over to the last spin. 

Note that a long distance singlet between the end spins 
may be realized in a ground state when those spins are very 
weakly coupled (J' = e/y/N « 1/\Z~N) to a spin chain 12 . 
This static approach to generate entanglement relies on cou- 
plings which are so weak that they can merely be regarded as 
perturbations. Such entanglement is not robust against ther- 
mal fluctuations due to the smallness of the gap (oc J' 2 = 
e 2 /N) between the ground state and a triplet state between 
the end spins. On the other hand our approach enables to 
generate entanglement dynamically even for J' as high as 
J'opt ~ l/(l°g^) 2 - Even when temperature is increased, the 
entanglement is not seriously disrupted till KbT exceeds the 
Kondo temperature (oc l/£ = 1/{N — 2)) after which the 
KC does not form. As a result, while in the dynamical ap- 
proach KbT < 1/(N — 2), in the static approach one has 
KbT < e 2 /N: thus, the long distance entanglement gener- 
ated through the dynamical approach is thermally more sta- 
ble. For instance, for e ~ 1CP 1 , our dynamical approach is 
robust for temperatures 100 times higher than those required 
for the static approach. In Fig. [2}l, we plot E m - as obtained 
in both approaches- vs. temperature for N = 10. In the static 
approach, the ground state is replaced by the thermal state, 
whereas in our dynamic approach it is the initial state which 
is taken to be the relevant thermal state. We ignore thermaliza- 
tion and relaxation during dynamics since the dynamical time 
scale, set by t opt , is fast enough (this is also an advantage over 
slow dynamical schemes with weak couplings^). 
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FIG. 3: (Color online) a) Different \GSi) for entanglement gen- 
eration through quench dynamics. Top: The ground states with 
£ < N/2 (no entanglement). Middle: N/2 < £ < N - 2 (some 
entanglement). Bottom: ^ — N — 2 (optimal entanglement), b) 
Decoupling the impurity from the chain. 




FIG. 4: (Color online) E m vs. iV after decoupling the first impurity. 



There are systems where the Kondo cloud mediated long 
distance entanglement may be observed such as spin chains in 
ion traps 2 -, with trapped electrons^, in chains of P donors in 



5»=a and Josephson chains with impuritie 
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III. CONCLUSIONS 



We have shown that substantial long range distance 
independent entanglement can be engineered by a non- 
perturbative quenching of a single bond in a Kondo spin chain. 
This is the first example where a minimal local action on a 
spin chain dynamically creates long range entanglement. In 
contrast to all known schemes for entanglement between in- 
dividual spins in spin chains, here the entanglement attains a 
constant value for long chains rather than decaying with dis- 
tance. We showed that, in the KR, the entanglement between 
the ending spins is mediated by the the KC. As the coupling 
is non-perturbative in strength (J' opt ~ l/(log-/V) 2 ), the en- 
tanglement generation is both fast and thermally robust. In 
the KR, the entanglement is periodic in time with a period of 
2t opt - this is a curious instance of exciting a long lived regu- 
lar oscillation in a gapless system. Both, the long distance en- 
tanglement mediated by the cloud and the periodic dynamics, 
provide a clear signature of the existence of the KC in a quan- 
tum system with impurities. These features are absent when 
the KC is absent, e.g., in the gapped DR. Beyond the spin 
chain based implementations, in principle, in double quantum 
dots one may detect the length of the KC by inducing and 
probing the entanglement oscillations between the dots. Our 
analysis evidences that interesting applications to quantum in- 
formation may arise exploiting relevant non-perturbative co- 
operative phenomena of condensed matter physics. 
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